Fundamental constants and tests of theory in Rydberg states of hydrogen-like ions 
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Comparison of precision frequency measurements to quantum electrodynamics (QED) predictions for Ryd- 
berg states of hydrogen-like ions can yield information on values of fundamental constants and test theory. With 
the results of a calculation of a key QED contribution reported here, the uncertainty in the theory of the energy 
levels is reduced to a level where such a comparison can yield an improved value of the Rydberg constant. 
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Quantum electrodynamics (QED) makes extremely accu- 
rate predictions despite the "mathematical inconsistencies and 
renormalized infinities swept under the rug" flj. With the as- 
sumption that the theory is correct, it is used to determine val- 
ues of the relevant fundamental constants by adjusting their 
values to give the best agreement with experiments J2). In 
this paper, we consider the possibility of making such compar- 
isons of theory and experiment for Rydberg states of cooled 
hydrogen-like ions using an optical frequency comb. We find 
that because of simplifications in the theory that occur for Ry- 
dberg states, together with the results of a calculation reported 
here, the uncertainty in the predictions of the energy levels 
is dominated by the uncertainty in the Rydberg constant, the 
electron-nucleus mass ratio, and the fine-structure constant. 
Apart from these sources of uncertainty, to the extent that the 
theory remains valid, the predictions for the energy levels ap- 
pear to have uncertainties as small as parts in 10 17 in the most 
favorable cases. 

The CODATA recommended value of the Rydberg con- 
stant has been obtained primarily by comparing theory and 
experiment for twenty-three transition frequencies or pairs of 
frequencies in hydrogen and deuterium [2]. The theoretical 
value for each transition is the product of the Rydberg con- 
stant and a calculated factor based on QED that also depends 
on other constants. While the most accurately measured tran- 
sition frequency in hydrogen (1S-2S) has a relative uncer- 
tainty of 1.4 x 1CP 14 [3], the recommended value of the Ryd- 
berg constant has a larger relative uncertainty of 6.6 x 1CP 12 
which is essentially the uncertainty in the theoretical factor. 
The main source is the uncertainty in the charge radius of the 
proton with additional uncertainty due to uncalculated or par- 
tially calculated higher-order terms in the QED corrections. 
This uncertainty could be reduced by a measurement of the 
proton radius in muonic hydrogen [4|, or by a sufficiently 
accurate measurement of a different transition in hydrogen. 
On the other hand, for Rydberg states, the fact that the wave 
function is small near the nucleus results in the finite nuclear 
size correction being completely negligible. Also, for Ryd- 
berg states, the higher-order terms in the QED corrections are 
relatively smaller than they are for S states, so theoretical ex- 
pressions with a given number of terms are more accurate. 

Circular Rydberg states of hydrogen in an 80 K atomic 
beam have been studied with high precision for transition 
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FIG. 1: Graph showing values of Z and approximate n that give a 
specified value of the frequency for transitions between states with 
principal quantum number n and n — 1 in a hydrogen-like ion with 
nuclear charge Z. Frequencies in the near infrared and visible range 
are indicated in color. 

wavelengths in the millimeter region, providing a determina- 
tion of the Rydberg constant with a relative uncertainty of 
2.1 x 10 -11 |6). With the advent of optical frequency 
combs [7 |, precision measurements of optical transitions be- 
tween Rydberg states have now become possible using fem- 
tosecond lasers. An illustration is the laser spectroscopy of 
antiprotonic helium |8|. Figure [T] gives iso-frequency curves 
corresponding to the spacing between adjacent Bohr energy 
levels (n to n — 1) in the 2-dimensional parameter space of 
the principal quantum number n and the nuclear charge Z 
for hydrogen-like ions. Much of this space is accessible to 
optical frequency synthesizers based on mode-locked fem- 
tosecond lasers, which readily provide ultra-precise reference 
rulers spanning the near-infrared and visible region of the 
optical spectrum (530 nm-2100 nm). Diverse techniques in 
spectroscopy (such as double-resonance methods) broaden the 
range of applications. Even when the absolute accuracy is lim- 
ited by the primary frequency standard (a few parts in 10 16 ), 
optical frequency combs can enable relative frequency mea- 
surements with uncertainties approaching 1 part in 10 19 over 
100 THz of bandwidth |9). The precisely controlled pulse 
train from a femtosecond laser can also be used directly to 
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probe the global atomic structure, thus integrating the optical, 
terahertz, and radio-frequency domains [ 10 1 . 

There are simplifications in the theory of energy levels of 
Rydberg states of hydrogen-like ions that, in some cases, al- 
low calculations to be made at levels of accuracy comparable 
to these breakthroughs in optical metrology. In the following, 
we write the known theoretical expressions for the energy lev- 
els of these ions, describe and give results of a calculation that 
eliminates the largest source of uncertainty, and list the largest 
remaining sources of uncertainty. We also make numerical 
predictions for a transition in two different ions as illustra- 
tions, look at the natural line width, and discuss what might 
be learned from comparison of theory and experiment. 

In a high-n Rydberg state of a hydrogen-like atom with nu- 
clear charge Z and angular momentum I = n — 1, the prob- 
ability of the electron being within a short distance r from 
the origin is of order (2Zr / nao) 2n+1 / (2n + 1)1, where do is 
the Bohr radius. Due to this strong damping near the origin, 
effects arising from interactions near or inside the nucleus are 
negligible, including the effect of the finite size of the nucleus. 

For I > 2, the theoretical energy levels can be accurately 
expressed as a sum of the Dirac energy with nuclear motion 
corrections _Edm> relativistic recoil corrections -Err, and ra- 
diative corrections .Bqed: E n = £ DM + -Err + Eqed- Re- 
views of the theory and references to original work are given 
in El ITTI [121 . The difference between the Dirac eigenvalue 
and the electron rest energy is proportional to 
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where a is the fine-structure constant, 5 = \k\ — 
^k 2 - (Za) 2 , k = (-l) l+3+1/2 (j + 1/2) is the Dirac spin- 
angular quantum number, and j is the total angular momen- 
tum quantum number. The energy level, taking into account 
the leading nuclear motion effects, but not including the elec- 
tron or nucleus rest energy, is given by lfT2ll 
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where h is the Planck constant, c is the speed of light, = 
a 2 m e c/2h is the Rydberg constant, tn = m c /m^ is the ratio 
of the electron mass to the nucleus mass, and fi r = 1 /(l + rjj) 
is the ratio of the reduced mass to the electron mass. 

Relativistic corrections to Eq. (|2]) associated with motion of 
the nucleus are classified as relativistic-recoil corrections. For 
the states with I > 2 under consideration here, we have 
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where lnfco(n, 1) is the Bethe logarithm. We assume that 
the uncertainty due to uncalculated higher-order terms is 
Za In (Za)~ 2 times the contribution of the last term in 
Eq. (|3J. 

Quantum electrodynamics (QED) corrections for high-Z 
states are summarized as 
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where a c is the electron magnetic moment anomaly and 
G(Za) is a function that contains higher-order QED correc- 
tions. Equation Q contains no explicit vacuum polarization 
contribution because of the damping of the wavefunction near 
the origin. Also in that equation, the uncertainties in the 
theory of a c may be eliminated by using the experimental 
value a c = 1.159 652 180 85(76) x 10~ 3 obtained with a one- 
electron quantum cyclotron lfl3l . 

The leading terms in G(Za) are expected to be of the form 
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The coefficients indicated by the letter A arise from the 
one-photon QED corrections; Aq$ and A%\ arise from the 
self energy, and Ago arises from both the self energy and the 
long-range component of the vacuum polarization. The term 
has been calculated for many states with I < 8, but not 
for higher-^ states before this work. The uncertainty intro- 
duced by this term if it were not calculated, based on plausi- 
ble extrapolations from lower-Z known values, would be the 
largest uncertainty in the theory and larger than the uncer- 
tainty from the Rydberg constant. The higher-order coeffi- 
cient A$i and the self-energy component of the coefficient 
A so are not known, but can be expected to be small. The 
vacuum polarization contribution to A 8 q is known 1 14 1 and is 
extremely small. The coefficient Bqo arises from two-photon 
diagrams and has not been calculated for high-/ states, but a 
comparison of calculated values of Beo [15 \ and A§o [ 16 1 for 
I < 5, suggests it has a magnitude of roughly 4A 60 , which 
is used as the associated uncertainty. Further, a term propor- 
tional to (a/n)BQi \n(Za)- 2 , that is nonzero for S and P 
states, vanishes for higher-/ states [17|. The term Cqq is ex- 
pected to be the next three-photon term, in analogy with the 
two-photon terms. 

In order to eliminate the main source of theoretical uncer- 
tainty in the energy levels, we have calculated the value of 
Ago for a number of Rydberg states. This calculation uses 
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TABLE I: Calculated values of the constant Aeo- The numbers in 
parentheses are standard uncertainties in the last figure. 
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methods from field theory, i.e., nonrelativistic QED (NRQED) 
effective operators which facilitate the calculation QjQ , and 
methods from atomic physics to handle the extensive angu- 
lar momentum algebra in the higher-order binding corrections 
of near-circular Rydberg states. Distinct contributions to the 
self-energy from high- and low-energy virtual photons, are 
matched using an intermediate cutoff parameter [19|. For 
near-circular Rydberg states, the radial wave functions have at 
most a few nodes, yet the calculation of A§q coefficients for 
these states is much more involved than for low-lying states. 
The reason is that in using the Sturmian decomposition of 
the hydrogen Coulomb Green function, as done for lower-n 
states, the radial integrations lead to sums over hypergeomet- 
ric functions with high indices, which in turn give rise to an 
excessive number of terms. For states with n = 8, there are 
of order 10 5 terms in intermediate steps, which is roughly two 
orders of magnitude more terms than for states with n = 2 
||20l . This trend continues as n increases making calculation 
at high n with this conventional method intractable. 

Here we report that the calculation has been done with a 
combined analytic and numerical approach based on lattice 
methods by using a formulation of the Schrodinger-Coulomb 
Green function on a numerical grid ETI . Provided quadru- 
ple precision (~ 32 significant digits) in the Fortran code is 
used, and provided a large enough box to represent the Ryd- 
berg states on the grid is used, the positive continuum of states 
can be accurately represented by a pseudo-spectrum of states 
with positive discrete energies. With this basis set, the virtual 
photon energy integration can be carried out analytically for 
each pseudo-state using Cauchy's theorem. This solves the 
problem of the calculation of the relativistic Bethe logarithms 
without the need for the subtraction of many pole terms, which 
would otherwise be necessary if the virtual photon energy 
were used as an explicit numerical integration variable. The 
results of this calculation for a number of states with n = 13 
to 16 are given in Table [I] 

We incorporate the results for A^o to numerically evaluate 
the theoretical prediction for the frequency of the transition 
between the state with n = 14, I = 13, j — lr and the 
state with n = 15, I = 14, j = ^ in the hydrogen-like ions 
He" 1 " and Ne 9+ . The constants used in the evaluation are the 
2006 COD ATA recommended values El . with the exception 
of the neon nucleus mass m( 20 Ne 10+ ) which is taken from the 



TABLE II: Transition frequencies between the highest-j states with 
77. = 14 and n = 15 in hydrogen-like helium and hydrogen-like 
neon. 



Term 


4 He+ i/(THz) 


20 Ne 9+ i/(THz) 


-Edm 


8.652 370 766 008(58) 


216.335 625 5746(14) 


-Err 


0.000 000 000 000 


0.000 000 0001 


-Bqed 


-0.000 000 001894 


-0.000 001 1841 


Total 


8.652 370 764114(58) 


216.335 624 3907(14) 



neon atomic mass [23], corrected for the mass of the electrons 
and their binding energies. Values of the various contributions 
and the total are given as frequencies in Table [II] Standard 
uncertainties are listed with the numbers where they are non- 
negligible. The theory is sufficiently accurate that the largest 
uncertainty arises from the Rydberg frequency ci?^, which 
is a factor in all of the contributions. There is no uncertainty 
from the Planck constant, since v = — E\±)/h. 

Table [TIT] gives sources and estimates of the various known 
uncertainties in the theory. To put them in perspective, in 
hydrogen, the relative uncertainty from the two-photon term 
£?6o for the 1S-2S transition is of the order of 10~ 12 due to 
disagreement between different calculations, whereas in the 
n = 14 to n = 15 Rydberg transition it is likely to be roughly 
5 x 10~ 19 , based on the smallness of the calculated value of 
the coefficient. The improved convergence of the expan- 
sion of the QED corrections in powers of Za is indicated by 
the fact that A^o is smaller by a factor of about 10 6 for the 
Rydberg states than the value Aqo « —30 for S states. 

The QED level shift given by Eq. Q is understood to be 
the real part of the radiative correction, while the complete 
radiative correction to the level £qed = -Eqed — i-T/2 is 
complex and includes an imaginary part proportional to the 
rate A = T/fr for spontaneous radiative decay of the level 
to all lower levels. For the highest-/ state with principal quan- 
tum number n, the dominant decay mode is an El decay to the 
highest-/ state with principal quantum number n — 1 ll24l . For- 
mulas in Ref. l24l give the nonrelativistic expression for the 
decay rate, which can also be derived from the nonrelativistic 
limit of the imaginary part of the level shift |25 1. 

As a first approximation, for transitions between states with 
quantum numbers n and n— 1 the ratio of the transition energy 



TABLE III: Sources and estimated relative standard uncertainties in 
the theoretical value of the transition frequency between the highest- 
j states with n — 14 and n = 15 in hydrogen-like helium and 
hydrogen-like neon. 
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Rydberg constant 
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FIG. 2: Graph showing values of Z and approximate n that give a 
specified ratio of the transition frequency to the natural width of the 
transition resonance between circular states with principal quantum 
number n and n — 1 in a hydrogen-like ion with nuclear charge Z. 



to the width of the line, is given by 
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where the expression on the right is the asymptotic form as 
n — > oo of the nonrelativistic value. Figure [2] gives a con- 
tour plot of the values of n and Z that give a specified value 
of Q based on the nonrelativistic result in Ref. 11241 . This is 
just a rough indication, since transitions with smaller I values 
will generally have a smaller Q, whereas transitions with a 
change of n greater than 1 will have a larger Q. The effect 
of possible asymmetries of the line shape on the apparent res- 
onance center has been shown to be small by Low 11261 . of 
order a(Za) 2 Eqeu- For the 1S-2S transition in hydrogen, 
such effects are indeed completely negligible at the current 
level of experimental accuracy ll27l . However, for Rydberg 
states of hydrogen-like ions, particularly at higher-Z, asym- 
metries in the line shape, some of which depend on details of 
the experiment, may be significant, and can be calculated if 
necessary. 

Recent advances in atomic-molecular-optical physics have 
generated an array of tools and techniques useful for engi- 
neering highly simplified atomic systems ||28l . In particular, 
observations of cold antihydrogen production at CERN illus- 
trate two ways for a cooled ion/antiproton to capture an elec- 
tron/positron in high-Z Rydberg states, either by three-body re- 
combination or by charge exchange [29 1. Properties of atomic 
cores have also been studied using a double-resonance detec- 
tion technique to observe the fine structure of Rydberg states 
produced by charge exchange in a fast beam of highly-charged 
ions ll30l . Using electron cooling l29l (and charge exchange), 
cold hydrogen-like ions can be recombined in high-^ Rydberg 
states from a variety of bare ions extracted from sources such 
as an electron beam ion source/trap (EBIS/T). Although two- 
photon spectroscopy is possible in certain cases, if the ions 
are confined in a trap within a region smaller than about half 



the wavelength of the radiation exciting the transition, Dicke 
narrowing also eliminates the first-order Doppler shift BP . 
Assuming T = 100 K, the relative second-order Doppler shift 
is about 3.5 x 10" 12 for He+ and 7 x 10" 13 for Ne 9+ . Tem- 
peratures in the range 4 K < T < 77 K are obtainable in 
cryogenic ion traps by resistive cooling [31] and by electron 
or positron cooling ll29l . For lower temperatures (T < 1 K), 
sympathetic laser cooling methods can be used PP . 

Of the variety of (n, I, Z) combinations of hydrogen-like 
ions, circular Rydberg states of low-Z ions seem the most fa- 
vorable for a comb-based determination of the Rydberg con- 
stant. On the other hand, some perturbations are smaller and 
line-widths are larger in heavier ions. Hence, using ions with 
a variety of (n, Z) combinations could be useful for exper- 
imental optimization and consistency checks, as well as for 
extending diversity of experiments used to determine funda- 
mental constants and test theory. 
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